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We present measurements of weak gravitational lensing cosmic shear two-point statistics using 
Dark Energy Survey Science Verification data. We demonstrate that our results are robust to the 
choice of shear measurement pipeline, either ngmix or imSshape, and robust to the choice of two- 
point statistic, including both real and Fourier-space statistics. Our results pass a suite of null 
tests including tests for B-mode contamination and direct tests for any dependence of the two- 
point functions on a set of 16 observing conditions and galaxy properties, such as seeing, airmass, 
galaxy color, galaxy magnitude, etc. We furthermore use a large suite of simulations to compute the 
covariance matrix of the cosmic shear measurements and assign statistical significance to our null 
tests. We find that our covariance matrix is consistent with the halo model prediction, indicating 
that it has the appropriate level of halo sample variance. We compare the same jackknife procedure 
applied to the data and the simulations in order to search for additional sources of noise not captured 
by the simulations. We find no statistically significant extra sources of noise in the data. The overall 
detection significance with tomography for our highest source density catalog is 9.7cr. Cosmological 
constraints from the measurements in this work are presented in a companion paper [T]. 


I. INTRODUCTION 

Cosmic shear, the weak gravitational lensing of galax¬ 
ies due to large-scale structure, is one of the most statisti¬ 
cally powerful probes of Dark Energy, massive neutrinos, 
and potential modifications to General Relativity nig. 
Due to its powerful potential as a cosmological probe, 
many ongoing and future surveys (Kilo-Degree Survey: 
KiD^ Hyper Suprime-Cam survey: HSCq the Dark 
Energy Survey: DEq^ the Large Synoptic Survey Tele¬ 
scope: LSSTQ Eucli(^and WFIRStQ will employ cos¬ 
mic shear as one of their principle cosmological probes. 
Cosmic shear two-point measurements, in their simplest 
form, are made by correlating the shapes of many mil¬ 
lions of galaxies as a function of their separation in angle. 
Additionally, if the galaxies can be separated as a func¬ 
tion of redshift, then tomographic cosmic shear measure¬ 
ments can be made by cross-correlating galaxies at differ¬ 
ent redshifts, which can probe the evolution of large-scale 
structure. The galaxies themselves have intrinsic shapes 
that are an order of magnitude larger than the cosmic 
shear signal, which means that cosmic shear measure¬ 
ments involve extracting small correlations from a large, 


* Corresponding author: beckermr@stanford.edu 
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shape noise-dominated background. Competitive cosmo¬ 
logical constraints from cosmic shear will require of order 
percent level or better measurements at all steps of the 
analysis, from shear measurement to the measurements 
of cosmic shear two-point functions (see, e.g., Weinberg 
et al. [4] or Kilbinger [5] for a review). 

Cosmic shear was first detected in 2000 [SHS]- The 
most recent results have detected correlated shapes on 
scales from a few to 60 arcminutes from the Deep Lens 
Survey [lO], the Sloan Digital Sky Survey HHHg, KiDs 
m and the Canada-France-Hawaii Legacy Survey mi, 
including in 6 redshift bins |15j . Future cosmic shear 
measurements will be very high signal-to-noise and over 
much larger survey areas, yielding a wealth of cosmolog¬ 
ical information. 

Cosmic shear measurements are challenging for a vari¬ 
ety of reasons. First and foremost, shear measurements 
are subject to biases that can arise from a number of 
sources. These biases are usually split into additive and 
multiplicative components. Sources of additive biases in¬ 
clude inaccuracies in the modeling of the point spread 
function (PSF), inaccuracies in correcting for the effect 
of the PSF on galaxy images, astrometric errors, and 
contaminating flux from nearby galaxies. Multiplicative 
biases can arise from the effects of noise on the shear 
measurement process, incorrect estimates of the size of 
the PSF, and, for model-fitting methods, mismatches 
between an object’s true underlying structure and the 
model employed in the shear measurement process. Ad¬ 
ditionally, many modern shear measurement methods re¬ 
quire accurate estimates of the distribution of galaxy 
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shapes and profiles in the absence of lensing to either 
serve as priors in the extraction of shapes from the data 
or to directly make corrections to the data. These priors 
can be estimated from high-resolution Hubble Space Tele¬ 
scope imaging, but must be matched to the observational 
sample under consideration. 

Significant computational and scientihc challenges in 
cosmic shear measurements remain, even in the pres¬ 
ence of perfect shear measurements. The cosmic shear 
field is the result of lensing by the non-linearly evolved 
matter density field. Accurate predictions for the non¬ 
linear matter power spectrum, even just for pure dark 
matter models, are computationally expensive and are 
needed at every point in parameter space in order to ex¬ 
tract cosmological parameters. Emulators, like the Coy¬ 
ote Universe [16] . have solved this problem for typical 
cosmologies and Dark Energy models, but neglect impor¬ 
tant physical effects, like galaxy formation, on the mat¬ 
ter power spectrum. Additionally, some physical effects 
of galaxy formation break the assumption that galaxies 
are randomly oriented in the absence of lensing. These 
effects, called intrinsic alignments, can introduce corre¬ 
lations in the shapes of galaxies that are not due to lens¬ 
ing, complicating the interpretation of cosmic shear mea¬ 
surements [see, e.g., dzi nH]. Furthermore, even if the 
mean signal can be modeled properly, the covariance ma¬ 
trix of cosmic shear measurements is dominated by sam¬ 
ple variance, requiring either extensive suites of numer¬ 
ical simulations or complicated halo model calculations. 
The (mis-)estimation of photometric redshifts (photo-zs) 
from imaging data is yet another important source of bias 
in the modeling of cosmic shear measurements. Finally, 
for precise cosmic shear measurements, lensing magnifi¬ 
cation, second-order lensing effects, and source selection 
effects will be important. 

In this work, we present cosmic shear measurements 
from Dark Energy Survey (DES) Science Verification 
(SV) data (Gruendl et al. in preparation; Rykoff et al. in 
preparation) using the shear catalogs by Jarvis et al. [T^ . 
We employ a combination of two shear estimation codes 
and two photometric redshift estimation codes, each of 
which takes a different approach to many of the issues de¬ 
scribed above. Additionally, we use a suite of ray-traced 
weak lensing simulations to compute the sample variance 
contributions to the covariance matrix of our measure¬ 
ments. We then present an extensive suite of tests of 
both the signals in the data and the covariance matri¬ 
ces. These tests include comparisons of the covariance 
matrices to halo model predictions, null tests of B-mode 
contamination, and null tests based on comparing the sig¬ 
nal between halves of the source galaxy sample split by 
survey metadata, like seeing, depth, etc. Overall, we find 
no statistically significant contamination. This paper is 
closely related to three other papers, namely the presen¬ 
tation of the DES SV shear catalog [19] , the presentation 
of the DES SV photometric redshifts for weak lensing 
[20] , and a companion paper that presents constraints on 
cosmological parameters using the measurements in this 


paper [T]. 

This work is organized as follows. In Section |ll| we de¬ 
scribe the DES SV shear catalogs and photometric red¬ 
shifts. Then we describe the mock catalogs used in this 


work in Section III Next, in Section |IV[ we present our 
detections of cosmic shear with DES SV data and our 
real-space two-point function estimators. Appendix 
describes alternate two-point estimators besides the real- 
space correlation functions used for the bulk of this work. 
We discuss the estimation and validation of our covari¬ 
ance estimation in SectionjVj Then, we describe our suite 
of null and consistency tests of our measurements in Sec¬ 


tion jVIj Finally, we conclude in Section VII The shear 
correlation functions and simulation covariance matrices 
from this work are available as online supplementary ma¬ 
terial with this paper. 


II. DATA 

The DES SV data with weak lensing measurements 
consists of 139 square degrees of hve-band imaging with 
roughly 7 exposures per band on average [2TH24] . The 
depth of the data is somewhat shallower than the ex¬ 
pected ~10-exposure average depth of the DES five-year 
data. The basic reductions and co-add source detection 
were done with the DES data management (DESDM) 
system as described in Desai et al. |2S] and Gruendl et 
al. (in preparation). We use the shear measurements 
from Jarvis et al. (THI performed on the DES SV Gold 
sample of galaxies (Rykoff et al., in preparation). For 
more information on the shear measurements and rec¬ 
ommended cuts, we refer the reader to Jarvis et al. [Ej. 
The shear measurement pipelines and photo-zs used in 
this work are described below for completeness. We use 
the “reduced shear” ellipticity definition [^. Finally, 
note that the two shear measurement codes used in this 
work are not identical, employing different cuts and dif¬ 
ferent parts of the DES SV data. Thus they have different 
overall source number densities and photometric redshift 
distributions. These differences, which we expect to be 
smaller in future DES analyses (see Jarvis et al. [ 19 ]), 
have no effect on the major conclusions of this work and 
are in fact important in verifying the robustness of our 
results. 


A. Shear Measurement Pipeline 1: ngmix 

The NGMi>([^ pipeline |17| uses sums of Gaussians to 
represent simple galaxy models [28] . The model param¬ 
eters of each object are sampled using Markov Chain 
Monte Carlo (MCMC) techniques applied to a full likeli¬ 
hood which forward models the galaxy and its convolu¬ 
tion with the PSF. The total likelihood for each object 


^ https://github.com/esheldon/ngmix 
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is a product of the likelihoods of the individual images 
of each object. The r-, i- and z-bands are ht simulta¬ 
neously with the same model shape, but different am¬ 
plitudes. The samples of the likelihood are then used 
with the lensfit algorithm to measure the shear 
of each object using a prior on the intrinsic distribution 
of shapes from the GREATS [ 30 ] release of the COSMOS 
galaxy sample. The final effective source number den¬ 
sity of the NGMix catalog is ~ 6.1 galaxies per square 
arcminute|^ Each source has an associated weight and 
we use the average sensitivity over both directions, as 
described in Jarvis et al. m 

B. Shear Measurement Pipeline 2: imSshape 

The iM3SHAPE[^pipeline is built on the imSshape code 
described in Zuntz et al. ED, with conhguration and 
modifications for its application to DES SV data 
described in Jarvis et al. [T9|. imSshape is a 
forward-modelling maximum likelihood code that uses a 
Levenberg-Marquardt algorithm to ht (in the conhgura¬ 
tion used here) two different models to galaxy images, 
one a de Vaucouleurs bulge and the other an exponen¬ 
tial disc, including the effect of the PSF and pixelization. 
The better-htting model is then used to give an ellipticity 
estimate. Maximum-likelihood parameter sets computed 
by imSshape and similar codes have a bias we refer to as 
noise bias ( 33 ] [ 33 ] . This bias is removed using a calibra¬ 
tion scheme based on the work of Kacprzak et al. |S3j . 
The scheme is applied to an ensemble of galaxies using 
the mean bias calibration for the ensemble; different sub¬ 
sets of objects thus use different correction factors. The 
hnal imSshape catalog has an effective number density 
of ~ 4.1 galaxies per square arcminute. Each source in 
imSshape has a weight, two additive noise bias correc¬ 
tions (one each for ei and 62 ) and a single multiplicative 
correction. 


C. Photometric Redshifts 

Based on an extensive comparison of four photo-z 
methods’ impacts on the two-point correlation function 
in Bonnett et al. [20j and a comparison of a much larger 
set of photo-z methods in Sanchez et al. |S4] . we have 
selected SkyNet [33 [S^ for our fiducial photo-z tomog¬ 
raphy. Galaxies are split into tomographic bins of equal 


® We use the following definitions of effective source density rieff 
and the effective shape noise per component ctsn? which are 
appropriate for the two-point function estimators employed in 
this work. rieff = = 

^ where Wi are the weights, Si 

are the sensitivities, a are the shear components, Q is the survey 
areas and the index i runs over all of the galaxies. 

® https://bitbucket.org/joezuntz/im3shape 


lensing weight for the NGMix catalog according to the 
mean of the photo-z PDF for each galaxy produced from 
SkyNet. The resulting tomographic bin boundaries are 
then used for galaxies in both shear catalogs. For a given 
shear code, the redshift distribution of each tomographic 
bin is estimated from summing the redshift probability 
distributions of each individual galaxy according to their 
weights assigned by the shear code. The relative agree¬ 
ment between the photo-z estimates and its impact on 
the correlation function is discussed in more detail in 
Bonnett et al. [20] . 


III. MOCK CATALOGS 

We use a set of 126 mock catalogs to compute the co- 
variance matrix of the shear correlation functions, E/B- 
mode statistics, power spectra and null statistics de¬ 
scribed in the following sections. These mock catalogs 
are constructed from seven sets of simulations consist¬ 
ing of three N-body light cones pieced together along 
the line of sight. We use 1050 /i“^Mpc, 2600 /i“^Mpc 
and 4000 h-^Mpc boxes with 1400^, 2048^ and 2048^ 
particles respectively. We use a flat, ACDM model with 
rim = 0.286, Da = 0.714, n, = 0.96, h = 0.7, Df, = 0.047, 
w = —1 and (Tg = 0.820. The initial conditions are gen¬ 
erated at redshift 49 with 2LPTic, a second-order La- 
grangian perturbation theory initial conditions generator 
[37] using linear power spectra from the GAME Boltzmann 
code [38]. The N-body evolution is computed with an 
efficient dark-matter-only version of the Gadget-2 code 
[39] . LGadget-2. We have implemented our own on-the- 
fly light cone generator directly into the LGadget-2 code 
(Busha et al. in preparation). We produce a full-sky 
light cone which formally replicates the N-body box eight 
times. However, each final simulation covers only one oc¬ 
tant of the full-sky, ~ 5, 000 square degrees, eliminating 
the replications. As the DES SV area with weak lensing 
measurements is only 139 square degrees, we divide each 
simulation into 18 different pieces using the observed SV 
mask to construct 126 total mock catalogs. This proce¬ 
dure has the advantage of properly computing the halo 
sample variance contributions to the lensing covariance 
matrices due to the fact that each patch is embedded in 
the large-scale modes of the box. 

We place lensing sources randomly in angle with in 
the DES SV mask (see Jarvis et al. US] for the details 
of the mask), and with the redshift distribution of the 
tomographic bins defined above. Then the weak lensing 
shear for each source is computed using the CALCLENS 
ray-tracing code m- In this application of CALCLENS, 
we use the pure spherical harmonic transform version 
with Nside=8192. Appendix [^ presents tests of the un¬ 
derlying simulations in comparison to simple expecta¬ 
tions from fitting functions to the matter power spec¬ 
trum. We find that the simple expectations from matter 
power spectrum fitting functions agree with the simula¬ 
tion to within sample variance, but that some resolution 
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FIG. 1. The measured shear correlation functions (,+/- for a single tomographic bin for the ngmix shape catalog (left) and 
imSshape shape catalog (right). The single tomographic bin corresponds to redshift distribution shown in Figure]^ z « 0.3 —1.3. 
Note that the redshift distributions of the two catalogs are not identical, so that the shear correlation functions are not expected 
to match. A detailed comparison of the two catalogs is described in Section [VI B[ Negative measurements are shown as upper 
limits. The error bars show the Ict uncertainties from the mock catalogs with the appropriate level of shape noise for each 
shear pipeline. The black solid lines show the predictions from a flat, ACDM model described in Section HI- not chosen to 
fit the data. 


issues remain on small scales. Note, however, that these 
small scales are excluded from the companion cosmologi¬ 
cal analysis [T] and that despite the resolution issues, we 
find excellent agreement between the covariances com¬ 
puted from the mock catalogs and the halo model, as 
discussed below. Thus for purposes of computing covari¬ 
ance matrices, the mock catalogs we have constructed 
are sufficient. Future work may require higher-resolution 
shear fields for covariance estimation. 


Finally, we generate the shape noise and other proper¬ 
ties in the mock by randomly drawing from the observa¬ 
tions separately for each tomographic bin. Importantly, 
we draw the intrinsic shape of each mock shear source 
separately from its other properties, like signal-to-noise, 
size, etc. Properties which have intrinsic spatial depen¬ 
dence in the survey (e.g. seeing, airmass, etc.) are drawn 
from the nearest real galaxy to each mock galaxy. See 
Section VIC for more details. These procedures ran¬ 
domise the shear field in the data and ensure that the 
mock catalogs have no correlations between the system¬ 
atic parameters and the shear field. 


IV. MEASUREMENTS OF COSMIC SHEAR 
TWO-POINT STATISTICS 


In this work, we focus on cosmic shear measurements 
made with two-point statistics, which are detailed in 
the following sections. A companion paper [T] presents 
the associated cosmological parameter constraints using 
these measurements, which use the real-space two-point 
correlation functions as the fiducial two-point estima¬ 
tor. We summarize results from alternate estimators in 
Section IV C and Appendix Note that although the 
choice of which two-point statistic to use is somewhat 
arbitrary, the companion cosmological analysis of this 
data [T] demonstrates that the exact choice of two-point 
statistic does not change the cosmological parameter con¬ 
straints from this data in a statistically significant way. 


A. Real-space Two-point Function Estimators 

We follow Miller et al. [4T] and estimate the two-point 
functions with 

= ^+±Wx (1) 

w,Wj{e - c)i,+/x(e - c)j-+/x 
'Yhi j WiWjSiSj 




X 


-t/x 
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FIG. 2. The measured shear correlation functions ^+/_ times Q in six angular bins and three tomographic bins for the ngmix 
shape catalog (left) and imSshape shape catalog (right). The tomographic bins correspond to those shown in Figure 
2 « 0.30 — 0.55,0.55 — 0.83,0.83 — 1.30, and are labeled from 1 to 3, increasing with redshift. Thus, panel ‘3-2’ shows the 
cross-correlation between the highest and middle redshift bins. The error bars show the Icr uncertainties from the mock catalogs 
with the appropriate level of shape noise for each shear pipeline. As in Figurethe black solid lines show the predictions from 
our hducial ACDM model — not chosen to fit the data. 


where z, j index the galaxies in the two sets we are cor¬ 
relating. Here e+/x are the estimated shears from the 
lensing analysis projected into the -I- (tangential) and x 
(cross) components rotated into the reference frame con¬ 
necting each each pair of galaxies {i,j} in the sum. The 
Wi are weights applied to each galaxy (typically inverse 
variance weighting; see Sec. [IT for each lensing code). 
The Si are multiplicative noise bias and/or lensing sen¬ 
sitivity corrections that are applied to the shears. We 
follow Miller et al. [41] and apply these corrections to the 
entire population of shears as opposed to applying them 
to each shear individually. We compared several differ¬ 
ent methods for incorporating the sensitivities into the 
two-point function estimator and find that they differ by 
at most ^ 2%. The are the additive bias corrections 
used for imSshape and are identically zero for NGMix 
per the definition of the lensfit method |29|. Finally, 
we use TreeCorip^ [42] to compute the shear correlation 
functions. 


https : //github. com/rinjarvis/TreeCorr 


B. Real-space Correlation Functions 


The real-space correlation functions without tomog¬ 
raphy are shown in Figure [^ We show ngmix on the 
left and imSshape on the right, with in the top rows 
and in the bottom rows. Negative measurements 
are shown as upper limits. The redshift distribution of 
sources for the non-tomographic analysis is shown in the 
top panel of Figure]^ for the SkyNet code. It appears in 
Figure[2that the correlation function may approach a 
constant value at large scales. Interestingly, Jarvis et al. 
[ig find that the mean shear across the survey for NGMix 
and imSshape is « 7 — 10 x 10“^. This level of mean 
shear would produce a constant floor in the shear correla¬ 
tion functions of « 5 —10 x 10“^. For the DES SV survey, 
the root-mean-square mean shear just due to shape noise 
and cosmic variance is « 4 x 10“'*. Thus it is not clear 
if this feature is an indication of systematic effects or a 
few sigma fluctuation in the mean shear due to a real 
physical effect. However, in the cosmological analysis of 
this data, all data points above 60 arcminutes were 
cut to avoid systematics in the PSF models [Il[ig. Thus 
we do not explore this issue further in this work. 


We generate estimates of the Icr uncertainties for each 
measurement by computing the covariance of the two- 
point functions over the simulation mock catalogs de¬ 
scribed in Section HI These mock catalogs are built sep- 
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FIG. 3. The estimated redshift distributions from SkyNet for the ngmix shape catalog (left) and the imSshape shape catalog 
(right). The full n(z) for objects with mean redshifts in the redshift range 0.3 < z < 1.3 (top) and the n(z) for three 
tomographic bins (bottom) are shown. The redshift distributions are estimated by summing and rescaling the photometric 
redshift probability distributions for each galaxy in the tomographic bin using the weights applied to the shear catalog. 


arately for each shear catalog in order to match the non- 
tomographic redshift distribution of the sources. The 
correction factor described in Hartlap et al. [13] is then 


applied to produce an unbiased estimate (see Section V A 
for a further discussion of the statistical properties of the 
covariance matrix estimate from the mock catalogs). The 
significance of the resulting measurement is then calcu¬ 
lated from this covariance as 


S/N = 


^dataC 


Ddel 


model^ 


^^model 


( 2 ) 


where C~^ is the inverse covariance matrix estimated 
from the mock catalogs, Idata is the vector of real-space 
shear two-point function measurements from the data, 
and ^modei is the vector of real-space shear two-point 
function measurements predicted from the cosmological 


model given above in Section III This quantity corre¬ 
sponds to the signal-to-noise of a least-squares estimate 
of a scaling parameter comparing our measurements to 
the theoretical model. This signal-to-noise measure will 
be an underestimate if the model employed is not well 
matched to the data. However, given the good match of 
our fiducial model to the data as shown in Figures and 
the degree to which the signal-to-noise is underesti¬ 
mated is small in this case. We use the COSMOSIS pack- 
ag^^by Zuntz et al. [U] to compute the shear correlation 
functions with the Takahashi et al. |45j non-linear power 
spectrum fitting function. See the companion paper [T] 
presenting cosmological constraints from these measure¬ 
ments for additional details on the model correlation 


https://bitbucket.org/joezuntz/cosmosis 


function ^modei computation. The covariance matrix has 
been validated through comparisons to both a detailed 
halo model prediction and jackknife estimates in single 
mock patches versus the survey data, which are discussed 
in detail in the Section]^ We find non-tomographic cos¬ 
mic shear detections at 6.5a and A.7a significance for 
NGMIX and im3shape respectively. 

Figure shows the full three-bin tomographic shear 
correlation function measurements for ngmix on the left 
and imSshape on the right. The redshift distributions 
of the three tomographic bins for the SkyNet code are 
given in the lower panels of Figure]^ In order to com¬ 
pute the covariance matrix of these measurements, we use 
the same procedure in the mock catalogs as for the non- 
tomographic case, except that we use the tomographic 
redshift distributions to assign the mock galaxies to dif¬ 
ferent tomographic bins. We additionally draw the shape 
noise in the mock from only the galaxies in the data in 
the same tomographic bin. We find overall tomographic 
cosmic shear detections of 9.7cr and 7.0cr for NGMIX and 
imSshape, respectively. Note that the ngmix catalog has 
more sources and extends to slightly higher redshift on 
average, yielding higher significance detections of cosmic 
shear. We have chosen three tomographic bins as a com¬ 
promise between gaining signal-to-noise in the data and 
having too many data points in order to use the mocks 
to compute the covariance matrix of the data. 

In Figures and the solid black line shows the ex¬ 
pected amplitude and shape of the shear correlation func¬ 
tions in the cosmological model given above. This curve 
is not a fit, and is merely presented as a reference for 
comparison. Due to the fact that the two catalogs have 
different redshift distributions, a direct comparison of the 
shear correlation functions between the two catalogs is 
not possible without further work matching the two cat- 















alogs and accounting for the shared shape noise, sample 
variance, and image noise between the two catalogs. This 
matched comparison is described further in Sec. |VIB| 


C. Alternative Two-point Statistics 


In Appendix we describe results from two alter¬ 
native two-point statistics of the shear field. These in¬ 
clude the methods of: (i) Becker and Rozo [35], which 
use a weighting of the real-space correlation estimates to 
construct efficient estimates of the Ci values and (ii) a 
second estimation of the spherical harmonic shear power 
spectrum using PolSpic Note that these es¬ 

timators weight the data at different angular scales dif¬ 
ferently than the default two-point correlation functions 
so that we do not expect to get identical results in terms 
of the significance of the cosmic shear detection. We do 
find detections of cosmic shear that are consistent with 
the conventional real-space estimators we use by default, 
indicating no strong preference for any given estimator. 
Tests of B-mode statistics from these estimators are dis¬ 
cussed in Sec. VIA where we again find consistency be¬ 
tween different two-point function estimation methods. 


V. ESTIMATING AND VALIDATING THE 
COVARIANCE MATRIX 

In this section, we present our covariance matrix and a 
set of validation tests. The fiducial covariance matrix for 
our measurements is estimated from the mock catalogs 
presented in Section |lll| First we compare the covariance 
matrix from the mock catalogs to halo model computa¬ 
tions. Second, we compare jackknife covariances in the 
data to the jackknife covariance computed from the mock 
catalogs. This procedure allows us to look for additional 
sources of noise and correlations in the data that are not 
present in the mock catalogs. 


A. Simulation and Halo Model Comparison 

We compare the covariance matrix computed from the 
simulations to that obtained from a halo model in Fig¬ 
ure m The simulation-based covariance matrix is com¬ 
puted by populating the mock catalogs with shear sources 
as described above in Section |TTl| and then computing the 
covariance of the measurements performed on the full en¬ 
semble of mock catalogs. The halo model covariance was 
computed with the CosmoLike covariance module (see 
Eifler et al. [49] and Krause et al. [50] for details). Fur¬ 
ther details of our halo model computation and the full 
tomographic covariance matrix are given in Appendix [C] 


http;//www2.iap.fr/users/hivon/software/PolSpice/ 


Briefly, we include the Gaussian, non-Gaussian and halo 
sample variance terms [e.g., El] and compute the halo 
model covariance at the same cosmology and with the 
same redshift distribution as was used in the mock cata¬ 
logs. 

We compare the general structure of the mock (up¬ 
per triangle) and halo model (lower triangle) covariance 
in the left panel of Figure [^ which shows part of the 
correlation matrix. Here we have shown a subset of the 
full set of tomographic bin combinations. The full corre¬ 
lation matrix is shown in Appendix [^ The right panel 
compares the amplitude of the two covariances by plot¬ 
ting the variance. Overall, we find good agreement in 
both structure and amplitude. 

We quantitatively test the agreement using a Fisher 
matrix computation. We compute the expected error 
on the degenerate parameter combination (78(11^/0.3)°'®, 
where as is the RMS amplitude of the linear matter 
power spectrum at redshift zero in a top hat window 
of 8 /i“^Mpc and flm is the matter density in units of 
the critical density at z = 0. This combination of pa¬ 
rameters is typically the best constrained by low-redshift 
cosmic shear data sets like the DES SV data. The exact 
degeneracy is computed in the companion cosmological 
constraints paper to this work [1]. We use the standard 
Fisher matrix formalism for cosmic shear [see, e.g., [55] 
and the same cosmological model as described above. We 
vary only the spectral index rig, as and flm in the Fisher 
matrix. 

We find that the error bars on cr8(^m/0.3)°'® from 
the halo model and mock covariances agree to approx¬ 
imately 10% without tomography, with the halo model 
yielding larger parameter uncertainties. When repeat¬ 
ing the same exercise with tomography, we find a larger, 
« 35% disagreement in the error bars, with the mocks 
yielding larger errors. However, we expect fluctuations 
in the uncertainties in parameters computed with the 
simulations due to the finite number of realizations used 
for the covariance computation. Dodelson and Schnei¬ 
der [53] estimate that this effect, in the Gaussian limit, 
increases the variance in the parameter estimates by a 
factor of 

{Nd-Np){N,-Nd-2) 

“ {N, - m - l){Ns - Nf, - A) 

where is the number of data points, Ng is the number 
of simulations and Np is the number of parameters. This 
factor is « 1 -I- Nd/Ns in the limit that ^ Nd ^ Np. 
Thus we expect a fractional uncertainty in the parameter 
uncertainties of « \/a — 1. In our case with tomography, 
Nd = 72, Ns = 126 and Np = 1. With these numbers, we 
get that the fractional uncertainty in the parameter un¬ 
certainty is « 118%. Thus the disagreement of « 35% we 
find with the halo model with tomography is not statis¬ 
tically significant. Without tomography, we find a frac¬ 
tional uncertainty in the uncertainty of « 56%, again 
indicating consistency. 

Importantly, these numbers are the fractional uncer¬ 
tainty in the uncertainty. For parameter estimates, the 
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FIG. 4. Comparison of the tomographic shear correlation function correlation matrix estimated from the mock catalogs and 
calculated from the halo model. The left plot shows the correlation matrix from the mock catalogs (upper left) and halo model 
(lower right). We show only the components for the first and last tomographic bins, plus their cross correlations. On the right, 
we show the square root of the diagonal elements of both covariance matrices, sorted in reverse numerical order. The open 
symbols show the results from the halo model and the closed symbols show the results from the mock catalogs. 



FIG. 5. Jackknife covariances in the mock catalogs and the data for the non-tomographic shear correlation functions. The 
left panel shows the correlation matrix using jackknifes in the data for ngmix (top left) and when averaging the jackknife 
covariances in the 126 mock catalogues (bottom right). The bottom left quadrant contains the correlations, the top right 
the correlations and the off diagonal components contain the cross-correlations. For each submatrix of the full correlation 
matrix, the angular scale increases from 2 arcminutes to 300 arcminutes. On the right, we show the diagonal elements of the 
jackknife covariance matrix in the data for ngmix (points) and when averaged over 126 mock catalogues (line). The grey band 
shows the standard deviation of diagonal elements over the 126 mock catalogs. 


fractional increase in the uncertainty on the parameter, 
equal to is the relevant quantity. For tomography, 
this fractional increase is « 55% and without tomography 
it is « 15%. Furthermore, we have assumed that the to¬ 
mographic analysis uses all 72 data points. As described 
in The Dark Energy Survey Collaboration tt al. [I], only 
36 of the 72 data points are used for tomography, bring¬ 
ing the fractional increase in the error due to the finite 
number of realizations down to only « 18%. Similar cuts 
are made for the non-tomographic analysis, using only 16 


of the 30 data points. This number of data points results 
in a fractional increase of the parameter uncertainties of 
only « 7% for the non-tomographic analysis. 

B. Jackknife Comparisons to Data 

While our mock catalogs include both sample variance 
and shape noise contributions, any spatially varying sys¬ 
tematic effects, like errors in the shear calibration, should 
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be included in the covariance matrix of the shear corre¬ 
lation functions as well. To search for these potential ef¬ 
fects, we use the jackknife covariance matrix of the shear 
correlation functions as a statistic to be compared be¬ 
tween the data and the mock catalogs. Any additional 
sources of noise in the data, which are captured by the 
spatial scale of our jackknife regions, will show up as a 
difference between the jackknife covariance as computed 
in the data versus the mock catalogs. 

We estimate the jackknife covariances from the data 
and our mock catalogs as follows. We divide both the 
mock catalog and data into 100 spatial sub-regions, em¬ 
ploying the k-means algorithm]^ These regions are then 
used to perform jackknife resampling. For the details of 
jackknife covariance estimation for cosmic shear correla¬ 
tion functions, we refer the reader to a (technical) com¬ 
panion paper where these choices are examined in further 
detail (Friedrich et al. [Mj, see also Norberg et al. for 
an application to galaxy clustering). We use the standard 
jackknife scheme, where all of the shear sources in an en¬ 
tire subregion are removed for each jackknife resampling, 
which is called the galaxy-jackknife in Friedrich et al. |54j . 

Note that we are not comparing jackknife covariances 
with the true covariances, but rather simply the co- 
variance in the shear correlation function across the DES 
SV survey to the same statistic computed with the mock 
catalogs. Thus the absolute correctness of the jackknife 
covariance matrix is not an issue for our test, since it is 
just a statistic that is sensitive to the effects for which we 
wish to search. The performance of empirical covariance 
measures for cosmic shear surveys is explored in Friedrich 
et al. [5^ . 

The comparison of our jackknife procedure between 
the mocks and the data is shown in Figure Here we 
plot the correlation matrix of the averaged jackknife co- 
variance from the 126 mock ngmix catalogs (left panel, 
on the bottom right) and the same computation in the 
DES SV data (on the top left). The right panel compares 
the diagonal elements of the jackknife covariance for 
and when averaged over 126 mock catalogs and when 
computed from the data for ngmix. Using the Fisher ma¬ 
trix procedure described above, we find that the error on 
CT8(r2m/0.3)°'® from the data jackknife covariance matrix 
agrees with the mean of the ensemble of errors on this 
parameter from the mock jackknife covariances to within 
one standard deviation of the error over the ensemble. 
Thus we conclude that there are no statistically signifi¬ 
cant sources of additional variance in the data compared 
to the mock catalogs. 


Implemented for python by Erin Sheldon, 
WWW.github.com/esheldon/kmeans_radec. 


VI. TESTS FOR RESIDUAL SYSTEMATIC 
ERRORS IN THE COSMIC SHEAR SIGNAL 

Systematic errors in shear measurements can be from 
a wide array of sources ranging from telescope optics and 
observing conditions to details in the modeling, measure¬ 
ment, and calibration of shapes. The development of 
tests to identify potential systematic errors is critical to 
verifying accurate measurement of cosmic shear. Toward 
this end, we devise a set of tests that should produce a 
null result when applied to true gravitational shear. The 
measurement of a significant non-zero result is then an 
indication of unresolved systematic errors in the shear 
catalog that could bias measurements. 

The DES SV shear catalogs have passed a rigorous 
set of both traditional and novel null tests that lay the 
groundwork for validating the precise measurements that 
will be made with ongoing DES measurements during 
the main survey observing period. These tests are per¬ 
formed both at the catalog level and during the process 
of validating specific measurements based on the shear 
catalogs. We describe the methodology and results of 
both traditional and new null tests for sources of poten¬ 
tial systematic errors in both the non-tomographic and 
tomographic measured cosmic shear signal in the next 
two sections. 

Catalog-level tests were performed by Jarvis et al. 
[B cf. their Section 8] and included tests for additive 
systematic errors related to spatial position, the PSF, 
and galaxy properties. These tests included the cross¬ 
correlation of the galaxies and the PSF. No significant 
additive systematic errors were found, and they put up¬ 
per limits on the potential additive systematic contribu¬ 
tion to in their Section 8.7. In addition, the over¬ 
all multiplicative bias of the shear estimates was tested 
with simulations in Jarvis et al. [19] . Jarvis et al. [19] 
concludes that both catalogs are consistent with having 
small overall multiplicative bias, but due to uncertain¬ 
ties in their ability to constrain this value, they suggest 
marginalizing over a prior on the multiplicative bias with 
a standard deviation of 0.05 (see Equation 8-12 of Jarvis 
et al. [IS])- This multiplicative systematic is treated in 
the cosmological analysis of this data [T], where it con¬ 
tributes to an increase in the uncertainties on the final 
cosmological parameters constrained with this data. 


A. B-mode Measurements 

The cosmic shear field can be characterized by E- and 
B-modes which differ in parity. At first-order in the grav¬ 
itational potential in General Relativity, cosmic shear 
produces a pure E-mode field [see, e.g., [SB]. However, 
contaminating signals, like that from the telescope point- 
spread function, tend to contain both E- and B-modes. 
Thus one of the first suggested tests of cosmic shear de¬ 
tections was verifying that the B-mode signal is consis- 
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FIG. 6. Tomographic B-modes in DES SV data for ngmix (left) and imSshape (right). The error bars are calculated from the 
simulation realizations using the shape noise appropriate for each catalog. The tomographic bins correspond to those shown in 
Fig. [3] and are labeled from 1 to 3, increasing with redshift. Thus, panel ’3-2’ shows the cross-correlation between the highest 
and middle redshift bins. The total x^/d.o.f., accounting for the correlations between the points in each panel, for ngmix is 
62.5/60 and for im3shape is 41.2/60. 


tent with zero [57] Many methods have been suggested 
for B-mode estimation [e.g.. [601165] . Here we use the esti¬ 
mators from Becker and Rozo |46] , which estimate band- 
powers using linear combinations of the shear two-point 
functions that optimally separate E- and B-modes [40] . 
These estimators are 



( 3 ) 


( 4 ) 


where the sums run over the angular bins of the shear 
two-point functions. The weight vectors /+/- are chosen 
to simultaneously minimize E- to B-mode mixing while 
also producing compact band-power estimates in Fourier- 
space. See Appendix [X] for more details. 

In Figure [^ we show a measurement of the tomo¬ 
graphic B-mode signal using the Becker and Rozo [15] 
band-powers. We find no statistically significant B- 
mode contamination, with a total y^/d.o.f. for NGMIX 
of 62.5/60 and for im3shape of 41.2/60. The error bars 
in this case are computed using the mock catalogs above. 
In Appendix]^ we verify this conclusion by computing a 


Small levels of B-modes are produced at second order in the 
gravitational potential, but these are small enough not to spoil 
the null test [see, e.g.. 15811591 . 


complementary measurement of the non-tomographic B- 
mode signal using an alternate estimation of the spher¬ 
ical harmonic shear power spectrum. We find the B- 
modes from this alternate technique are consistent with 
zero with a y^/d.o.f. = 4.5/7 for NGMix and 6.3/7 for 
im3shape. Finally, note that Becker and Rozo [15] band- 
power measurements of the non-tomographic B-mode sig¬ 
nal are presented in Jarvis et al. uni using the methods 
and mock catalogs of this work. The non-tomographic 
B-mode measurements were again found to be consis¬ 
tent with zero, with y^/d.o.f. = 22.3/20 for NGMix and 
16.1/20 for im3shape. 


B. Consistency Between the Shear Pipelines 

We further test for consistency between the shear cata¬ 
logs split into tomographic bins by selecting only sources 
which pass the selection cuts for both codes. For this 
subset of sources, we then compare the shear auto- and 
cross-correlation functions for each bin. Due to the fact 
that the two catalogs have the same sample variance, 
have similar shape noise and have correlated shear mea¬ 
surement errors, the error bars on the difference between 
the two correlation functions is much smaller than that 
on the correlation functions themselves. We account for 
this effect by constructing mock catalogs where a given 
mock galaxy is assigned its shape noise for each shear 
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FIG. 7. Difference over error in the tomographic correla¬ 
tion functions for matched shear catalogues from ngmix and 
imSshape. We show imSshape minus ngmix. The total 
X^/d.o.f. accounting for all correlations is 46.8/72. 


measurement code, ngmix or imSshape, from the same 
real galaxy. 

This comparison is shown in Figurej^for the shear cor¬ 
relation function for imSshape minus NGMixlHI We find 
that the shear correlation functions from the codes are 
statistically consistent over the full range of scales from 2 
to 300 arcminutes, giving a x^/d.o.f. = 46.8/72. Finally, 
note that this test is similar to the differenced shear cor¬ 
relation function test presented in Section 8.6 of Jarvis 
et al. [El- For their test, they examine the shear cor¬ 
relation function of the the difference in the NGMIX and 
imSshape shear estimates using the matched catalogs. 
They find that below «3 arcminutes, the catalogs do not 
meet the requirements for additive systematic errors, set 
by the expected precision of the cosmological constraints. 
The test presented in this work is generally less sensitive, 
but complementary, to the differenced shear correlation 
function. 


C. Two-Point Null Tests 

Even with a carefully chosen set of null tests at the 
catalog level, it is still possible that systematic errors. 


We have completed this test for the ratio of the shear correlation 
functions and without tomography, finding similar results. 


which can be due to complex interplays between differ¬ 
ent aspects of data and analysis, may influence the cos¬ 
mic shear measurement. To test for any uncorrected sys¬ 
tematic errors remaining in the measured cosmic shear 
signal, we attempt to measure the variation in as a 
function of survey and galaxy properties that may be 
correlated with sources systematic errors. For each sur¬ 
vey or galaxy property, the shear data is split in half, 
and the correlation functions of each half are compared. 
We use a reweighting method to ensure that the redshift 
distribution of each half is the same in order to remove 
any cosmological dependence from this null test. If the 
photo-zs and shear measurements are correct, then the 
shear correlation functions of the two halves should be 
consistent to within the noise of the shear measurements 
and the redshift reweighting. If they are not, this would 
indicate either uncorrected systematics, selection effects 
from the split, or non-shear differences in the two halves 
such as intrinsic alignments. 

Due to the fact that each half is drawn from the same 
area in many cases, the standard error bars computed 
for the shear correlation functions are not correct for 
this test. We instead use the mock catalogs described 
above to compute the error on the difference between 
the two halves relative to the full sample, accounting for 
shared sample variance, as described below. It is impor¬ 
tant to note that this is a simultaneous test of both the 
photometric redshifts and the shear calibrations. This 
feature is in fact desirable because both of these quan¬ 
tities can contribute to biases in the shear correlation 
functions. We have used both the survey property maps 
described by Leistedt et al. [66] and also properties di¬ 
rectly produced by the shape measurement codes. The 
16 various systematic parameters are described in Ta¬ 
ble |lj Finally, Jarvis et al. m found that making cuts 
on signal-to-noise and size could lead to a selection bias 
in the population of shear values due to preferentially se¬ 
lecting galaxies that look more or less like the PSF. We 
attempt to minimize this problem by using the “round” 
measures of signal-to-noise, size and surface brightness. 


1. Methodology 

The galaxies in each half-sample must be reweighted 
so that the total n(z), computed from summing the in¬ 
dividual p[z) for each galaxy according to its weight, 
matches between the two half-samples. Matching the 
redshift distributions of the two halves removes any cos¬ 
mological dependence in each null test. For the data, the 
extra weights are computed using Ridge Regression (or 
Tikhonov regularization) [67]. We use the Ridge Regres¬ 
sion algorithm to solve for an additional weight for each 
galaxy, which when used with the shear measurement 
weights described in Section|n]to compute the n(z), pro¬ 
duces a matching redshift distribution between the two 
half-samples. The Ridge Regression algorithm solves the 
linear least-squares problem with an additional regular- 
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FIG. 8. An example of the redshift re-weighting procedure used when comparing the correlation function between galaxies 
split into bins of galaxy or survey properties. Left: The SkyNet redshift distribution for each half of the ngmix data, split 
into upper (blue) and lower (red) bins of signal-to-noise ratio {S/N)r before (dashed) and after (solid) re-weighting, compared 
to the full sample n{z) (black solid curve). Right: The distribution of weights applied to each galaxy to produce the solid n{z) 
lines, generated as described in Sec. |VI C| 


ization parameter a, minimizing 

||Rv-t||-h ||a(v-I)|| (5) 

where ||...|| denotes the least-squares norm, R is the ma¬ 
trix of galaxy p{z) ’s each weighted by the leasing weights 
given in Section [III 
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for n galaxies and m photo-z bins with leasing weights Wi 
and galaxy p(zys pji, t is the target photo-z distribution, 
V is the vector of new weights for which we are solving 
and I is the identity vector. The parameter a governs 
the flexibility of the weight selection - the smaller the 
value, the better matched the reweighed n{z) are — and 
is adjusted to prevent a significant contribution of nega¬ 
tive or large weight values, which may impact the validity 
of the null tests. We find that a = 5 x 10“^^ produces 
an optimal match between the two half-samples while 
keeping the weights v sufficiently regular for our photo- 
zs and lensing weights. This value may not generalize 
to other lensing weights or photo-zs. We match the red¬ 
shift distribution of each half-sample to that of the full 
sample (i.e. t is the redshift distribution of the full sam¬ 
ple). This procedure is more stable than matching one 
half to another since smaller weights are needed for each 
half. The application of the Ridge Regression algorithm 
then produces a new weight v, which is combined mul- 
tiplicatively with the lensing weight in the calculation 


of the correlation functions. The resulting reweighting 
for galaxies split into bins of low and high galaxy detec¬ 
tion signal-to-noise for NGMix is shown in Fig. The 
left panel shows the n(z) for each half before (dashed) 
and after (solid) reweighting compared to the full sam¬ 
ple. The corresponding weight histograms are shown in 
the right panel. 


We use the 126 DES SV-shaped mock catalogs de¬ 
scribed above to compute the variance and significance of 
the differences between the shear correlation functions in 
each half-sample. In the mock catalogs, we select a sub¬ 
set of galaxies in narrow redshift slices to match the n(z) 
distribution for the full galaxy catalog. Random shape 
noise is generated from the shear catalog and applied to 
the mock catalogs, and the property with which we split 
the galaxy sample in half is then mapped onto the galax¬ 
ies in each mock via a nearest neighbour algorithm in 
angular position, and redshift. This preserves the same 
spatial patterns as exist in the data, but the shears have 
been randomised so that there is no correlation with this 
property. We then apply the same procedure to each 
mock as applied to the data to directly compute the er¬ 
ror bars on the difference via Monte Carlo, with the ex¬ 
ception of using the true mock point redshift values to 
reweight the n(z) histograms of each half instead of a 
p{z) estimate for each galaxy. We expect this difference 
will only underestimate the variance. Any statistically 
significant deviations then indicate that the there may 
be a residual systematic error in the shear catalogs re¬ 
lated to the quantity split upon, which has affected the 
measured two-point correlation function. 
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FIG. 9. Null tests for the ngmix two point correlation function based on a variety of catalog and survey properties as described 
in Table 1. Each panel for a given property shows the difference between the relative to its error for the galaxies in the 
upper and lower halves of the sample split into bins by the magnitude of the quantity. The two halves of the sample have been 
reweighted to have the same redshift distribution. The error on the difference is computed directly via the mock catalogs. Grey 
bands are shown representing the Icr and 2a variance at each value of 0. Adjacent points in angle are correlated. 


2. Results 

The split null tests on are presented in Figure for 
NGMIX and Figure for imSshape. This is repeated in 
Appendix [P] for For each quantity (panel), the differ¬ 

ence in is shown at each value of 9 relative to the Icr er¬ 
ror in the difference from the mock catalogs. Grey bands 
corresponding to la and 2 cr errors are shown for com¬ 
parison. The corresponding statistical significance of the 
null tests for imSshape and ngmix are given in Table 
We find that for both NGMix and imSshape the null tests 
pass with deviations smaller than 2cr (y^/d.o.f. = 17.8/8) 
for all tests except for NGMix airmass. Note that be¬ 
cause NGMIX has a higher source density, it is gener¬ 
ally more sensitive to residual systematic errors in these 
tests. While this detection is still weak, it warrants eval¬ 
uating whether this difference in the galaxy population 
halves will have a significant bias on the correlation func¬ 
tion. To test this, we also show in Table |l] the difference 
= ^+(upper) — C-i-(iower) relative to the Icr error on the 
full sample measurement. For NGMIX airmass, this differ¬ 


ence is approximately one-third of the statistical error on 
the measurements and consistent with the level of bias 
in several other quantities. Of slightly lesser significance 
are splits in the magnitude of NGMix PSF ei and 62 , for 
which PSF 62 has the largest difference in between 
upper and lower halves — though still small compared 
to the statistical error. 

There is some subtlety in interpreting the significance 
of these null tests. First, due to physical effects not ac¬ 
counted for in the simulations, some tests could yield 
non-zero results but not indicate systematic errors in the 
data analysis itself. For example, if the level of intrin¬ 
sic alignments differs between galaxies split by colour, 
then these null tests could fail and yet the shear mea¬ 
surements themselves could be free of systematic errors. 
Second, these tests could also, in principle, flag differ¬ 
ences between the shear calibrations of galaxies of dif¬ 
ferent types, which although interesting, may not ul¬ 
timately impact cosmological constraints from the full 
sample, which could be unbiased on average. Third, as 
stated above, it is not clear from these tests alone if any 
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FIG. 10. Null tests for the imSshape two point correlation function based on a variety of catalog and survey properties as 
described in Table 1. See Fig. [^for details. 


deviations are due to the shear measurements or the pho¬ 
tometric redshifts. Finally, note that the values from 
these tests are not independent, due to correlations in the 
underlying quantities used to construct the tests (e.g., the 
survey depth is correlated with the seeing). We have per¬ 
formed a large number of null tests, so to the extent that 
the values between many of the tests should be inde¬ 
pendent, we do expect some apparent deviations purely 
from statistical fluctuations. However, we have not at¬ 
tempted to combine the tests in order to quote an overall 
significance. 

VII. CONCLUSIONS 

In this work, we present cosmic shear two-point mea¬ 
surements from Dark Energy Survey Science Verification 
data. We find an overall detection significance of 9.7cr for 
our higher source density catalog, NGMix. We addition¬ 
ally present multiple advances in band-power estimation, 
covariance estimation, simulations versus theory, and null 
tests for shear two-point correlations. Through this work 
we demonstrate that our measurements are robust and 
free of statistically significant systematic errors. 


We demonstrate that the covariance matrices derived 
from the DES SV mock simulations presented in this 
work are consistent with the halo model, including the 
halo sample variance terms. We also compare the vari¬ 
ance in the mock catalogs to the variance in the DES SV 
data by comparing jackknife covariances computed in the 
data and mock catalogs. The structure of the covariance 
matrices is very similar and we detect no statistically sig¬ 
nificant sources of additional variance in the data. 

We find that the B-mode signals in the data are con¬ 
sistent with zero and that the two shear estimation codes 
agree well. We additionally present a set of simultane¬ 
ous null tests of the photo-zs and shear measurements, 
performed by splitting the shear sample in half according 
to some parameter and comparing the shear correlation 
functions of the halves. We find that these tests pass with 
no statistically significant indications of biases. We ex¬ 
pect null tests similar to those developed here to have in¬ 
creased utility in future cosmic shear analyses, where the 
statistical power is larger and the requirements for con¬ 
trolling systematic errors and shear selection effects are 
more stringent. The DES itself will have nearly « 36x 
more data and will measure cosmic shear at significantly 
higher signal-to-noise, so that these tests will be very 
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Property 

X" [d.o.f. = 8] 
NGMIX (imSshape) 

AC-l/<r(e+) 

NGMIX (imSshape) 

Description 

Signal-to-Noise 

4.9 ( 5.2 ) 

0.05 ( 0.49 ) 

Signal-to-noise of galaxy detection 

Galaxy Size 

5.3 ( 10.7 ) 

-0.3 ( 0.15 ) 

Galaxy size (deconvolved with PSF) 

Galaxy Colour 

7.3 ( 2.2 ) 

-0.31 ( -0.32 ) 

g — z colonr 

Surface Brightness 

7.8 ( 8.7 ) 

0.33 ( -0.32 ) 

Galaxy surface brightness 

RA 

7.0 ( 8.8 ) 

0.24 ( 0.28 ) 

Galaxy right ascension 

Dec 

4.0 ( 6.2 ) 

-0.24 ( -0.57 ) 

Galaxy declination 

E(B-V) 

5.1 ( 6.2 ) 

0.23 ( 0.06 ) 

Mean extinction 

Air Mass 

20.7 ( 13.8 ) 

0.31 ( 0.46 ) 

Mean r-band air mass 

Exposure Time 

4.7 ( 6.8 ) 

0.18 ( 0.3 ) 

Mean total r-band exposure time 

Mag. Limit 

4.4 ( 7.4 ) 

0.18 ( 0.45 ) 

Mean r-band limiting magnitude 

Sky Sigma 

1.7 ( 13.0 ) 

-0.02 ( -0.08 ) 

Mean r-band RMS sky brightness 

Sky Brightness 

5.0 ( 14.3 ) 

-0.05 ( -0.27 ) 

Mean r-band sky brightness 

FWHM 

6.4 ( 3.3 ) 

-0.23 ( -0.13 ) 

Mean r-band PSF FWHM 

PSF ei 

16.8 ( 13.5 ) 

0.12 ( -0.37 ) 

Galaxy PSF ei 

PSF 62 

17.1 ( 7.5 ) 

-0.58 ( -0.22 ) 

Galaxy PSF 62 

PSF Size 

2.6 ( 5.6 ) 

-0.1 ( 0.42 ) 

Galaxy PSF size 


TABLE I. Summary of null tests for ngmix and imSshape. Results are given as ngmix (imSshape). The values are given 
for the differences between the two-point correlation function calculated from galaxies that fall within one of two bins in each 
catalog or survey property. Also shown is the magnitude of the difference relative to the la error of the measurement of on 
the full sample. 


useful. 
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Appendix A: Alternative E- and B-mode Statistics 

In this appendix we consider alternative statistics of 
the shear held, verifying that our conclusions above, es¬ 


pecially that the B-modes are consistent with zero, do 
not depend on the choice of statistic. These alterna¬ 
tive statistics include the band-powers of Becker and 
Rozo [15] and power spectra band-powers estimated with 
PolSpic^^ [171 148j . 

1. Band-powers 

The band-powers of Becker and Rozo [15] use the 
methods of Becker [ID] to estimate Fourier-space band- 
powers directly from linear combinations of the real-space 
two-point functions. The hnal band-power estimates can 
be computed from the underlying E-mode power spec¬ 
trum as 

7 din £ £2 

E = J ^—CEE{i)W+{e) (Al) 

where IF-|-(£) is the band-power window function com¬ 
puted from the coefficients {f+i, f-i} in Eqs. I&EI 
See Becker and Rozo [15] for more details. The optimal 
computation of the band-powers requires computing the 
effective radial bin window functions of the shear corre¬ 
lation function points. Instead in this work we just use 
the geometric approximation to the bin window func¬ 
tions to compute the amplitudes {/+i,/_i}. This pro¬ 
cedure means that the band-powers do not separate E- 
and B-modes as well as they could in principle. However, 
when comparing to a fiducial cosmological model below, 
we do compute the band-power window function using 
estimates of radial bin window functions from the data. 
These window functions are computed via interpolating 
the weighted counts in each radial bin of the estimated 
shear two-point function. We have compared the results 
of this procedure for computing the window functions to 
estimates of the window functions from counts in finer 
bins. We find unsurprisingly that the bin window func¬ 
tions are quite smooth and thus the interpolation is ac¬ 
curate enough for our purpose. 


2. Spherical Harmonic Power Spectrum 

The cosmic shear power spectrum can also be esti¬ 
mated in spherical harmonic space, which has the advan¬ 
tage of being faster and less memory intensive than work¬ 
ing in real-space. In view of upcoming wide field galaxy 
leasing surveys, e.g. the full five year DES dataset, we 
therefore investigate the applicability of standard spher¬ 
ical harmonic space methods to weak leasing. For this 
purpose, we use the PolSpice [17113H] code together with 
the HEALPix [55] package, which has been applied to, 
amongst other things, CMB polarization data [e.g.. I69j. 


http;//www2.lap.fr/users/hivon/software/PolSpice/ 
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FIG. 11. Band-powers in DES SV data for ngmix (top) and 
imSshape (bottom). The error bars indicated by the grey 
bands are calculated from the simulation realizations using 
the shape noise appropriate for each catalog. The dotted lines 
show the band-power window functions W+ {tj scaled so that 
their peak values are 2 x 10“®. The solid line is the prediction 
for the shear power spectrum for the fiat, ACDM model given 
above. The dashed line shows the integral of the band-power 
window functions over the shear power spectrum. 



FIG. 12. Spherical harmonic shear power spectrum estimated 
using PolSpice. The left and right panels correspond to the 
NGMIX and imSshape catalogs, respectively. The top and bot¬ 
tom panels show the E- and B-modes, respectively. The mea¬ 
surement uncertainties are estimated using the mock catalogs. 
The black solid lines show the predictions for the flat, AGDM 
model given above. Note that the theoretical prediction has 
been convolved with the PolSpice kernels, which relate the 
true to measured power spectra. The S/N values for the E- 


modes are computed as outlined in Section 


IV A 


and the x 


values for the B-modes indicate consistency with zero. The 
reported values take into account correlations between the 
band-powers. 


PolSpice is based on the fast correlation function ap¬ 
proach described in [47] and [48] . The method is designed 
to exploit the advantages of both real and spherical har¬ 
monic space: to limit computation time and resources, 
the data are analyzed in spherical harmonic space. In or¬ 
der to facilitate demasking, the power spectrum is trans¬ 
formed to real-space in an intermediate step. In real- 
space, the survey mask can simply be corrected for, since 
the masked correlation function is the product of the un¬ 
masked correlation function and the correlation function 
of the mask. More precisely, the algorithm first calculates 
pseudo-C^’s from pixelized and masked galaxy ellipticity 
maps which are then transformed to the correlation func¬ 
tion. The real-space correlation function is then divided 
by the correlation function of the mask to correct for 
finite survey effects and inverted to obtain the full-sky 
power spectrum, removing E- to B-mode leakage in the 
mean. Incomplete sky coverage implies that the inver¬ 
sion can only be performed on angular scales for which 
the correlation function can be estimated thus introduc¬ 
ing Fourier ringing in the inversion process, which can be 
reduced by apodizing the correlation function. Both the 
apodization and finite integration range introduce kernels 
which relate the power spectra measured by PolSpice to 
the underlying true power spectra. These kernels can be 
computed for a given apodization scheme and integration 
range and can therefore be corrected for when comparing 
measurement to theory (for details see Chon et al. [48] ). 

For our analysis, we pixelize the galaxy ellipticities 
onto a HEALPix pixelization of the sphere with a resolu¬ 
tion of Nside=1024, where each pixel covers a solid angle 
of «11.8 arcmin^. In order to obtain a robust estimate of 
the shear field, we need to correct for multiplicative bias 
in the measured ellipticities. Since the correction factors 
described in Sections |II Aj and jll B| are noisy estimates of 
the true corrections, we determine the mean sensitivity or 
multiplicative bias correction for our galaxy sample and 
apply this mean correction to the pixelized maps. As the 
power spectrum is estimated from maps constructed from 
the discrete values of the galaxy ellipticities, we apply a 
conservative masking scheme to maximize galaxy num¬ 
ber density. We therefore adopt the DES SV LSS mask 
used for galaxy clustering measurements EOj. This mask 
is identical to the DES SV mask used for weak leasing 
except that it restricts analyses to the largest contiguous 
region overlapping the SPT-E held by selecting the area 
with 60 < ra [deg] < 95 and —60 < dec [deg] < —40. It 
further considers only regions with survey limiting mag¬ 
nitude in the i-band > 22.5 (i.e. all regions considered to 
provide at least lOcr measurements for objects at i-band 
= 22.5; [7D]). For the power spectrum measurement, we 
limit all integrations to scales smaller than 9max = 15 
degrees and we apodize the correlation function with a 
Gaussian window of 0fwhm = 10 degrees. Finally, we 
correct the measured power spectra for the HEALPix pixel 
window function and compress them into 7 band-powers 
with PolSpice band-power kernels. 

The noise power spectrum needs to be computed from 
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simulations. In order to produce noise only maps from 
the data, we remove correlations in the ellipticity maps 
by rotating each galaxy ellipticity by a random angle. 
We then estimate the noise power spectrum as the mean 
of the power spectra of 100 such random realizations. 
This procedure yields shape noise estimates consistent 

with Ci^sN = where (t^ pj^ is the variance of either 

component of the mean ellipticity per pixel and Upix is the 
angular number density of HEALPix pixels. Comparing 
the measured shape noise to the galaxy-based Gaussian 

shape noise estimate C^^sn = ’ where is the 

variance of either component of the galaxy ellipticities 
and Ugai denotes the galaxy number density, we find that 
the latter underestimates the measured shape noise. This 
suggests that the galaxy ellipticity distribution is non- 
Gaussian and the Gaussian approximation can therefore 
only be applied after averaging the galaxy ellipticities 
over pixels. We test the pipeline using Gaussian field 
realizations and the mock catalogs. 


3. Results 

Figure shows the non-tomographic band-powers us¬ 
ing the methods of Becker and Rozo [15], their window 
functions as the dotted lines, and their error bars com¬ 
puted with the mock catalogs as the grey bands. We 
find a detection significance G.lcr and 5.7cr for NGMIX and 
imSshape, respectively. These detection significances are 
similar to the real-space two-point functions. Finally, 
the solid line shows the expected shear power spectrum 
amplitude assuming the flat, AGDM model given above. 
The dashed line shows for each band-power the integral 
of the band-power window function over the shear power 
spectrum. 

Figure shows the results for the PolSpice statis¬ 
tics. We find a detection of cosmic shear of 5.6 (t and 5Aa 
for NGMIX and imSshape respectively for the PolSpice 
statistics. Note that the PolSpice statistics do not use 
as many high-£ modes as the real-space band-powers 
or the real-space correlation functions, so that one ex¬ 
pects a lower detection significance. We also find that 
the B-modes are statistically consistent with zero for the 
PolSpice statistics. 

Finally, note that these two estimators process the 


data in different ways (e.g., averaging the data in pixels 
versus computing real-space correlation functions), have 
different sensitivities to shot noise, and have different 
Fourier-space window functions. We thus do not expect 
them to give precisely the same results in Fourier-space 
for the shear power spectrum. However, we do expect 
that when treated self-consistently they should give 
statistically consistent results for cosmological parame¬ 
ters, as demonstrated in the accompanying cosmological 
analysis of this data [T|. 

Appendix B: Validation of the Mock Catalogs 

In this section we present a validation test for the mock 
catalogs. We compare the shear correlation functions 
measured in the mock catalogs in tomographic bins with 
the theoretical expectation from the Takahashi et al. jlSj 
fitting function for the matter power spectrum. The re¬ 
sult of this test is shown in Figure [T^ We find that 
at high redshift the small-scale shear correlation func¬ 
tions are suppressed relative to the theoretical expecta¬ 
tion. Note however that this numerical effect is below the 
scales where the two-point functions are being used for 
cosmological parameter estimation (« 2 — 4 arcminutes 
for and « 25 — 55 arcminutes for see Table 2 of 
The Dark Energy Survey Collaboration et al. [I])- Addi¬ 
tionally, we only estimate the covariance of the two-point 
functions from the mock catalogs, not the mean signal. 
Within the noise of our mock covariance matrix, the over¬ 
all parameter uncertainties are consistent when using the 
halo model versus the simulation covariance (see Sec. 
for a quantitative comparison). This fact may indicate 
that the covariance is less sensitive to these numerical 
effects than the mean signal. Future work may require 
higher-resolution shear fields for covariance estimation. 

Appendix C: Detailed Covariance Matrix Validation 

In this section, we present further details of the val¬ 
idation of the covariance matrices, including our tomo¬ 
graphic halo model computations and the comparison to 
the simulations. The halo model covariance was com¬ 
puted with the CosmoLike covariance module (see Eifler 
et al. [ini and Krause et al. [sni for details). 

In the halo model, the covariance of tomographic shear 
power spectra is given by |7TH75] 


Cov{C:^{h),C^J{h)) = 






r dH 

f \ 




2n3 
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with n* the number of source galaxies in tomography bin i, the ellipticity dispersion, A(li) = d^l ~ 
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FIG. 13. The shear correlation functions in the mock catalogs compared to the expected values from Takahashi et al. [45] for all 
three tomographic bins (labeled in the top left corner from left to right). In the top panel, the solid lines show the theoretical 
expectation, the bands show the Ict sample variance estimate and the dashed line shows the mean from the mock catalogs. 
is in red and is in blue. In the bottom panels, we show the fractional deviation of the mean signal in the mock catalogs from 
the expected values from Takahashi et al. [3^ in units of the sample variance. data below ~ 2 — 4 arcminutes and data 
below « 25 — 55 arcminutes is not used for the final cosmological analysis in The Dark Energy Survey Collaboration et al. [1] 
due to the expected baryonic effects in the matter power spectrum. 


2TTliAli the integration area associated with a power 
spectrum bin centered at h and width Ali, and 

and T^^usv convergence trispectrum of source red- 
shift bins i,j, k and I in the absence of finite volume ef¬ 
fects and the halo sample variance contribution to the 
trispectrum [^Eg. Our halo model implementation for 
these terms is described in Eifler et al. m- 

Note that Equation j CI j ignores the so called finite-area 
effect (cf. Sato et al. |75| or Friedrich et al. [54]), linear 
beat-coupling terms [e.g.,|75] and linear dilation terms Li 
et al. [e.g.,|77]. For a survey of the size of DES-SV the 
finite-area effect is expected to be negligible. Further¬ 
more, ignoring this effect is at most conservative since 


Cov(4^(0i),^^'(02)) = J ^Uo/4{l0l) J fJrJo/4 


it will slightly overestimate the statistical uncertainties. 
The beat-coupling terms are negligible compared to the 
halo sample variance terms (and even the non-Gaussian 
terms, see e.g., Takada and Jain [75]). Further, the lin¬ 
ear dilation terms reduce the effect of the beat-coupling 
terms and are negligible m- Finally, we have ignored the 
effects of masking (except for the total area of the survey 
in the halo sample variance terms). We have found with 
Gaussian simulations that the effects of the details of the 
mask, besides the overall survey area, are negligible when 
computing cosmological constraints. 

The covariance of angular shear correlation functions 
is then given by 


:'d2) Gov(Q^(;i),c^-'(i2)) (C2) 


where we use the results of Joachimi et al. |78| to simplify 
the calculation of the Gaussian part of the covariance. 


Figure shows the full tomographic correlation ma¬ 
trix, comparing the halo model on the lower-right and 
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FIG. 14. Comparison of the shear correlation function correlation matrix estimated from mock catalogs and calculated from 
the halo model. Figurej^shows a subset, those for tomographic bin combinations (1,1), (1,3) and (3,3), of the covariance matrix 
elements shown in this figure. The correlation matrix from mock catalogs is on the upper-left and that from the halo model is 
on the lower-right. 


the mock catalogs on the upper-left. The overall struc¬ 
ture of the covariance matrices is similar in both compu¬ 
tations, but the mock catalogs exhibit more noise in the 
off-diagonal components. 


Appendix D: Additional Two-Point Null Tests of 


We have repeated an identical analysis for to that 
described in Sec. |VI C| for We show the results of 


the tests for imSshape in Fig. flSjand for NGMix in Fig. 
[T6| Qualitatively, comparing to Figs. |^ & IT there is 
an indication that some of the larger deviations in the 
figures for may be due to additive systematic errors. 
For example, there is an offset in the difference of 
based on values of airmass at the 2cr level that disappears 
for The corresponding and difference values are 
given in Table [IT] There are no significant indications of 
systematic errors in these null tests for though this 
may simply be due to the poorer constraining power of 
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FIG. 15. Null tests for the ngmix two point correlation function based on a variety of catalog and survey properties as described 
in Table 1. See Fig. [^for details. 
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Property 

X^[d.o.f. = 8] A^_/ct( 5 _) 

NGMIX (im3SHAPE) NGMIX (im3SHAPE) 

Description 

Signal-to-Noise 

5.8 ( 1.8 ) 

-0.07 ( 0.03 ) 

Signal-to-noise of galaxy detection 

Galaxy Size 

2.5 ( 5.0 ) 

-0.23 ( -0.35 ) 

Galaxy size (deconvolved with PSF) 

Galaxy Colour 

7.1 ( 3.8 ) 

-0.3 ( 0.04 ) 

g — z colour 

Surface Brightness 

4.4 ( 5.2 ) 

-0.04 ( -0.06 ) 

Galaxy surface brightness 

RA 

2.9 ( 3.0 ) 

0.06 ( - 0.22 ) 

Galaxy right ascension 

Dec 

4.9 ( 3.5 ) 

-0.35 ( -0.37 ) 

Galaxy declination 

E(B-V) 

2.8 ( 4.9 ) 

- 0.22 ( - 0.02 ) 

Mean extinction 

Air Mass 

2.7 ( 3.4 ) 

- 0.01 ( -0.08 ) 

Mean r-band air mass 

Exposure Time 

4.5 ( 2.5 ) 

-0.35 ( 0.0 ) 

Mean total r-band exposure time 

Mag. Limit 

2.2 ( 3.3 ) 

-0.29 ( -0.43 ) 

Mean r-band limiting magnitude 

Sky Sigma 

3.8 ( 5.6 ) 

-0.21 ( -0.3 ) 

Mean r-band RMS sky brightness 

Sky Brightness 

4.0 ( 6.2 ) 

-0.27 ( -0.42 ) 

Mean r-band sky brightness 

FWHM 

4.1 ( 4.5 ) 

- 0.2 ( -0.08 ) 

Mean r-band PSF FWHM 

PSF ei 

2.7 ( 7.9 ) 

-0.37 ( -0.55 ) 

Galaxy PSF ei 

PSF 62 

6.8 ( 5.8 ) 

-0.5 ( -0.33 ) 

Galaxy PSF 62 

PSF Size 

1.2 ( 3.8 ) 

-0.08 ( - 0.1 ) 

Galaxy PSF size 


TABLE II. Summary of null tests for ngmix and imSshape. Results for ngmix and imSshape are given as ngmix (imSshape). 
The values are given for the differences between the two-point correlation function calculated from galaxies that fall within 
one of two bins in each catalog or survey property. Also shown is the magnitude of the difference relative to the Icr error of 
the measurement of on the full sample. 
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